Domain Wall Solutions of Spinor Bose-Einstein Condensates in an Optical Lattice 
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We studied the static and dynamic domain wall solutions of spinor Bose-Einstein condensates 
trapped in an optical lattice. The single and double domain wall solutions are constructed ana- 
lytically. Our results show that the magnetic field and light-induced dipolar interactions play an 
important role for both the formation of different domain walls and the adjusting of domain wall 
width and velocity. Moreover, these interactions can drive the motion of domain wall of Bose ferro- 
magnet systems similar to that driven by the external magnetic field or the spin-polarized current 
in fermion ferromagnet. 
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During the past several decades, the dynamics of spa- 
tial domain wall have attracted more attention in ferro- 
magnetism [1, 2]. The conventional ferromagnet is usu- 
ally composed of fermions which contribute the main site- 
to-site exchange interaction owing to the direct Coulomb 
interaction among electrons and Pauli exclusion princi- 
ple. This exchange interaction causes spin wave unstable, 
and its developing instability brings about the appear- 
ance of spatial domain walls or magnetic soliton [3| . The 
contribution of the magnetic dipole-dipole interaction to 
domain wall formation is usually neglected in practice be- 
cause it is typically several orders of magnitude weaker 
than the exchange coupling in fermion ferromagnet. 

However, the description of magnetism composed of 
Bosons are not well explored. Fortunately, a typical 
Bose system was realized in cold spinor ^^Rb gases [J, |5| 
and ^^Na gases [6] which provided a totally new envi- 
ronment to understand magnetism comprehensively. The 
ferromagnetic spinor Bose gases have attracted numerous 
theoretical and experimental interests [7-16], and some 
general properties as in conventional fermion ferromag- 
nets have been observed, such as spontaneous symmetry- 
broken ground state [8], spin domains and textures 
and normal spin-wave excitations spectrum Q. Espe- 
cially, the spontaneous symmetry breaking in ^^Rb spinor 
condensates ^1] clearly shows ferromagnetic domains 
and domain walls in a Bose ferromagnet. Recently, the 
ferromagnetism in a Fermi gas is also of interest [17] in 
the field of ultracold atoms. 

Spinor Bose-Einstein condensates (BECs) trapped in 
an optical potential have internal degrees of freedom due 
to the hyperfine spin of the atoms which liberate a rich 
variety of phenomena. When the potential valley is so 
deep that the individual sites are mutually independent, 
spinor BECs at each lattice site behave like spin mag- 
nets and can interact with each other through both the 
magnetic and the light-induced dipole-dipole interactions 
which are no longer negligible due to the large number of 
atoms at each lattice site. These site-to-site dipolar inter- 



actions can cause the ferromagnetic phase transition [9| 
leading to a "macroscopic" magnetization of the conden- 
sate array and the spin-wave like excitation [9, 10] and 
magnetic soliton [l8|,|l9| analogous to the ferromagnetism 
in soHd-state physics, but occur with bosons instead of 
fermions. Also, spinor BECs in an optical lattice is a 
typical physical realization of Salerno model [20]. 

In this paper, we explore the domain wall solutions of 
spinor BECs trapped in an optical lattice and the roles of 
the magnetic and light-induced dipole-dipole interactions 
for the type of domain wall solutions and domain wall 
width and velocity. 

We consider F = 1 spinor condensates trapped in a 
one-dimensional optical lattice formed by two-polarized 
laser beams counter propagating along the ?/-axis and the 
two laser beams are detuned far from atomic resonance. 
Under the tight-binding approximation, the Hamiltonian 
takes the form M-flOl 



(1) 



where is the nth collective spin operator, defined as 
Sn = al^{n)Fav0^v{i^)^ with a being annihilation operator 
and F being the vector operator for the hyperfine spin 
of an atom. The first term in the Hamiltonian results 
from the spin-dependent interatomic collisions at a 
given lattice site, with A'^ = (l/2)Aa / (i^r |0^(r)|^, 
where is proportional to the difference between 
the s-wave scattering lengths in the triplet and singlet 
channels [8], and (/>^(r) is the ground-state wave function 
for nth site. The direction of the magnetic field B 
is along z-axis and 7^ = fi^dF is the gyromagnetic 
ratio, with /i^ being the Bohr magneton and qf 
the Lande factor. The last two terms describe the 
site-to-site spin coupling induced by both the static 
magnetic and the light-induced dipole-dipole interac- 



tions [lo|, with the form J^' 
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-fUo/{24Ah^kl) JJ dvdv'fc iv') e^-^^-r^- . 
cos {kLV - kLy')cos{kLy)e^i • W(r') • e_i |0^(r - r')|^ 
|0^(r)|^ + 3J!^i/2. Here the cutoff function 

fc (r) = describes the effective interaction 

range of the hght- induced dipole-dipole interaction, 
with Lc being the coherence length associated with 
different decoherence mechanisms [21] and A being the 
detuned frequency. The wave number /cl = 27r/AL, the 
transverse coordinate = + is the width 

of the lattice laser beams, and /7o denotes the depth of 
optical lattice potential. The e±i are unit vectors in the 
spherical harmonic basis, and the form of tensor W (r) 
can be found in Ref. jTo|- Eq. ([T]) shows that the static 
magnetic and the light-induced dipole-dipole interaction 
can lead to the isotropic spin coupling denoted by J^®° 
and anisotropic spin coupling in the transverse direction 
denoted by J^^^^ • 

From Hamiltonian ([1]) we can derive the Heisen- 
berg^ equation of motion at /cth site for the spin, i.e., 
ihdSk/dt = [Sk^ H]. When the optical lattice is infinitely 
long and in the ultra-low temperatures for condensa- 
tion, the operator can be treated as a classical vector, 
Sk S{y,t). Here we assume all nearest-neighbor in- 
teractions are same, which is a good approximation in 
one-dimensional optical lattice j22|. Then we get the ef- 
fective Landau-Lifshitz equation [18| 



ls-^-Sx\2Jd'(^ 
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S'e,)+jBe,], (2) 



where J = 2J*''^" + J'"° and do = Al/2 being the lattice 
constant. In a rotating frame around 2:-axis with angular 
frequency ■j^B/h the spin vector S is related with the 
original one by the transformation S+ = S"^ + iS^ = 
gi^^-i-isBt/h^ Thus Eq. (0) becomes 
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where the superscript {') is omitted for pithiness, and 
the time t and coordinate y is in unit to = h/ (2J) and 
do, respectively. The parameter /3 = 4J*^^"/J can be 
controlled by tuning the lattice laser beams and the shape 
of the condensate at each lattice site. Such a character 
makes the lattice atomic spin chain an ideal tool to study 
a diversity of spin-related phenomena. 

Static domain wall solutions: We firstly seek for 
the single domain wall solution in the form 



cos 77q 

cosh ©0 



sm r]Q 
cosh 6o 



,5^ = tanhGo, (4) 



where 9o = koy — uoot + rj'o with t^q, t^q being two real 
parameters, and /cq, ujo to be determined. Substituting 
Eq. (H)) into Eq. ([3j) we find that cjq = and ko = a/— ^, 






FIG. 1: (Color online) The illustration of the double domain 
solutions (a) 5*^, (b) 5*^, and (c) in Eq. (|5|) with condition 
P = -0.2 and ki = 0.3. The two components S"" and 
possess of the double peak located at sinhGi = ibl/^i. 



and the solution (jH in fact denotes the static domain 
wall. This result implies the parameter /3 < which can 
be realized by a blue-detuned lattice (A < 0), where the 
condensed atoms are trapped at the standing- wave nodes 
and the laser intensity is approximately zero. As a result, 
the light-induced dipole-dipole interactions is very small 
and the spin coupling is governed mainly by the static 
magnetic dipole-dipole interaction which implies << 
1. In this case the domain wall width, defined by I/Zcq, 
can be adjusted mainly by J^®°. Under the condition 
jtran/ jiso ^ _o.25, the domaiu wah width is about 3do 
with do being the lattice constant. 

Eq. ([3|) has a norm form of Landau-Lifshitz type for a 
spin chain with an anisotropy. It can be solved by inverse 
scattering method where a couple of Lax equations is 
introduced for constructing the analytical solutions. In 
terms of our earlier results [23] we get two static double 
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FIG. 2: (Color online) The illustration of the double domain 
solutions (a) 5'^, (b) , and (c) in Eq. (|6|) with condition 
/3 = —0.2 and cp = 0.5. The two components and S'^ 
possess of the single peak located at z = 0. 



domain wall solutions. The first has the form 

= -2(5i(cos<l>isinh6i)/Ai, 
Sy = -2(5i(sin^isinhei)/Ai, 

= 1-2/Ai, (5) 

where Ai = 1 + ^?sinh^ei, Gi = [P/{Ski) - 2ki]y, 
$1 = [2ki^ /3/{Ski)]\ Si = {Wkf^/3) /{Wkf-/3), 
and ki is a real parameter. The illustration of domain 
wall solution in Eq. (|5]) is shown in Fig. 1. From Fig. 
1 and Eq. ([5]) we see that the two components and 
precess around the component with the frequency 
{2ki -\- /3/ and the z-component of spin vector 

is conservative. The component 5'^(or S^) possess of the 
double peak located at sinhOi = which varies 

with time periodically as shown in Fig. 1. The width 
of double domain wall is 1/ {/3/ (Ski) — 2ki), which is in- 
verse proportion to the parameter /3, while the maximum 
absolute value of valley for is constant. 



The other solution can be written as 

5^ = 2[cot sin ((^-c&2) cosh e;]/A;, 
sy = 2[cot cos ((/:)-<i>2) cosh e;]/A;, 
s' = 1-2/a;, 



(6) 



where A^ = 1 + cot^ cosh^ , ^2 = P [cos'^ (p) t, 
©i = — (sin ^) y , and (/? is a real parameter. The il- 
lustration of domain wall solution in Eq. (|6]) is shown 
in Fig. 2. Different from the former one, the two 
components and possess of the single peak lo- 
cated at = which also oscillates with time period- 
ically as shown in Fig. 2. The width of double do- 
main wall is l/(v'—/^ sin (/:)), which is inverse proportion 
to \J The maximum absolute value of valley for 
is 1 — 2/ (1 + cot^ (/:?) which is also independent on the 
parameter /3. We can rewrite 5*^ in Eqs. (|5|) and (|6|) as 
= {{Sl + 1) tanh^ 61 - l)/{{Sl - 1) tanh^ (-61) - 1) 
and (cot^ (/:)—! +tanh^ 02)/(cot^ tanh^ O'^), respec- 

tively, which implies the double domain wall solutions can 
be written asymptotically as a nonlinear combination of 
two single domain wall solution in Eq. (jH). 

For observation of the above domain wall, one of the 
possible physical realizations of a gas of dipolar BECs 
can be provided by electrically polarized gases of polar 
molecules or by applying a high dc electric field to atoms 
[2^ . In order to induce the dipole moment of the order 
0.1 D (Debye) one needs an electric field of the order of 
10 V/cm and the corresponding s-wave scattering length 
10-1000 A. 

Dynamic domain wall solutions: For red-detuned 
lattices, the condensed atoms are trapped at the max- 
ima of the intensity of the standing wave laser and the 
spin coupling is dominantly determined by the light- 
induced dipole-dipole interaction. In particular, the spin 
coupling is anisotropic in this case. By controlling the 
laser parameters, we may always make the light-induced 
dipole-dipole interaction dominate over the static mag- 
netic dipole-dipole interaction, i.e., J*^^" >> J^^^ > 
and P ^ 2. In this case we assume that the single mov- 
ing domain wall solution of Eq. (jSj admits the form 



tanhe2,S'^ 



cos r]i 
cosh 02 ' 



sin r]'{ 
cosh02 ' 



(7) 



where 62 = k2y + uj2t + r][, with r][^ rj'l being two real 
parameters, and /c2, ^2 to be determined. Substituting 
Eq. dll) into Eq. ^ we obtain 002 = P / 2 sin {2r]'{) and 
k2 = ±^/P smr]^{. The solution ([7]) shows that the dy- 
namics of domain wall is restricted in (/c2, co'2)-space, with 
I ^2 1 max = 1^2 max I = /^/2. The width of domain 

wall, defined by 1/ |/c2 1 , is inverse proportion to the square 
root of parameter /3 determined mainly by light-induced 
dipole-dipole interaction and the angle r]i of the three 
components of pseudospin vector S. The domain wall 
velocity, i.e., v = uj2/k2 = iy^cos?]'/, is dependent on 
the parameters /3 and r]^. When 7^'/ = nvr, n = 0, 1, 2, 
the domain wall velocity attains its maximum value 
Vmax = V?. When r]'{ = (n + l/2)7r,n = 0,1,2,..., the 
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FIG. 3: (Color online) The scattering process of two moving 
domain walls (a) , (b) , and (c) 5*^ expressed in Eq. 
(|8|), in which two domain walls propagate with the opposite 
direction and the same absolute value of domain wall velocity 
in the condition /3 = 1.5 and k, = tt/S. 



solution in Eq. ([7]) represents the static domain wall solu- 
tion of Eq. (|3|) similar the case of /3 < 0. It is interesting 
to estimate the velocity of domain wall. As a example, 
we consider the F = 1 electronic ground state of ^^Rb 
which the Lande factor gp = —1/2 and the gyromagnetic 
ratio 7^ = — As from Ref. [9] we estimate J^®° 
is about 1.1 X 10~^^J with = 852 nm. Asumming 
that J*^^" = 10J^^° we have P ^ 1.9 and the maximum 
domain wall velocity is about v = a/3<^o/^o ^ 130 nm/s. 
For chromium atoms, it has a magnetic moment and 
the corresponding domain wall velocity is about tens of 
domain wall velocity of ^^Rb with the same assumption. 

The above results show that the domain wall velocity 
can be controlled by the parameter P resulting from the 
external magnetic and light field-induced dipole-dipole 
interaction and the direction of pseudospin vector S. In 
fermion ferromagnet, a magnetic domain wall is a spa- 
tially localized configuration of magnetization, in which 



the direction of magnetic moments gradually inverses. 
When a spin-polarized electric current flows through a 
domain wah, the spin-polarization of conduction elec- 
trons can transfer spin momentum to the local magne- 
tization, thereby applying a spin-transfer torque, which 
manipulates the motion of domain wall [25|, |26| similar 
to that driven by applied an external magnetic field [13 • 
It is interesting that our domain wall solution in Eq. ([7]) 
shows the same properties as that in fermion ferromag- 
net, i.e., for red-detuned lattices the light-induced dipolar 
interactions can be seen as the external force to drive the 
motion of domain wall, and it has the potential contribu- 
tion for the research of domain wall motion in the Bose 
ferromagnet systems. 

The dynamic double domain wall solution can be con- 
structed by employing Hirota method [28] which is an 
effective straightforward technique to solve the nonlinear 
equations. Here we only mention the main idea of Hirota 
method briefiy. Firstly, it applies a direct transforma- 
tion to the nonlinear equation. Then, by means of some 
skillful bilinear operators the nonlinear equation can be 
decoupled into a series of equations. With some reason- 
able assumptions the exact solutions can be constructed 
effectively. Performing the normal procedure in Ref. ^] 
we get dynamic double domain wall solution as 

5^ = [e^cosh^(A:^^)-sinh^(4^)-cos^/>:]/A2, 
Sy = 2e^/^ cos cosh (^2^) cosh (cj2t)/A2, 

= -2e^/^ sin cosh (fc2^)sinh(cj2t)/A2, (8) 

where A2 = cosh^ (^2^) + sinh^ (^2^) + cos^ with the 
real parameter a^, /c2 = sin k,^ uj'2 — /3/2 sin (2/^), 

and 6> = In (cos^ n) . In figure 3 we plot the evolution 
of double domain wall solution S in Eq. (|8]). From 
figure 3 we see that the double domain solution in Eq. 
dHJ presents a general scattering process of two moving 
domain walls which propagate with the opposite direc- 
tion and the same absolute value of domain wall velocity 
y/PmsK. Analysis reveals that there is no amplitude ex- 
change among the components for two domain walls. 
However, there is a small shift for the center of each do- 
main wall during collision. 

In conclusion, we have presented two kinds of domain 
wah solutions of dipolar spinor BECs in an optical lattice 
based on an effective Hamiltonian of anisotropic pseu- 
dospin chain. In real experiment, the magnetic field and 
light-induced dipolar interactions are related with the de- 
tuned term which can be controlled easily by the external 
magnetic and light field. Our results show that we can 
control this site-to-site spin coupling of the condensates 
at each lattice to form the different domain wall solu- 
tions. Especially, since these magnetic field and light- 
induced dipolar interactions are highly controllable, the 
spinor BECs in an optical lattice as an exceedingly clean 
system offer a very useful tool to study spin dynamics 
in periodic structures and to understand ferromagnetism 
comprehensively. 
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